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Abstract 

A diabatic (configuration-fixed) constrained approach to calculate the potential energy surface 
(PES) of the nucleus is developed in the relativistic mean field model. As an example, the poten¬ 
tial energy surfaces of obtained from both adiabatic and diabatic constrained approaches 

are investigated and compared. It is shown that the diabatic constrained approach enables one 
to decompose the segmented PES obtained in usual adiabatic approaches into separate parts 
uniquely characterized by different configurations, to follow the evolution of single-particle orbits 
till very deformed region, and to obtain several well defined deformed excited states which can 
hardly be expected from the adiabatic PES’s. 
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I. INTRODUCTION 


The relativistic mean field (RMF) theory is one of the most successful microscopic 
models in nuclear physics. The RMF theory incorporates from the beginning very impor¬ 
tant relativistic effects, such as the existence of two types of potentials (Lorentz scalar and 
four-vector) and the resulting strong spin-orbit interaction, a new saturation mechanism 
by the relativistic quenching of the attractive scalar field, and the existence of anti-particle 
solutions. The Lorentz covariance and special relativity make the RMF tl^qry more ap- 

2, 3|. It has 
3|, exotic 


Q 


pealing for studies of high density and high temperature nuclear matter 
achieved success in describing many nuclear phenomena for stable nuclei 


nuclei 


as well as supernova and neutron stars 6|. The RMF theory provides a new 


explanation for the identical bands in superdeformed nuclei ^ and the neutron halo ^li 
predicts a new phenomenon — giant neutron halos in heavy nuclei close to the neutron 


drip line 


gives naturally the spin-orbit potential, the origin of the pseudospin sym¬ 


metry |9 |,Im as a relativistic symmetry 


11 


12 


13l |. and spin symmetry in the anti-nucleon 


spectrum Jl4, and well describes the magnetic rotation [l5|, the collective multipole ex¬ 
citations [l6| , and the properties of hypernuclei, etc. Lately, the ground-state properties 
of over 7000 nuclei have been calculated in the RMF-I-BCS model and good agreements 


with existing experimental data are obtained 


17l |. Recent and more complete reviews of 


the applications of the RMF model, particularly those on exotic nuclei, can be found in 

Refs. ybj. 

In order to describe the shape of the atomic nucleus and understand the fusion and 
fission processes, it is crucial to obtain the potential energy surface (PES) of the nucleus 


as a function of the deformation [20|. In phenomenological methods, the PES is obtained 
by minimizing the total energy of the system with some shape parameters. In micro¬ 
scopic models, such as Skryme Hartree Fock and RMF model, PES can be obtained in a 
constrained calculation. There are two different ways to obtain the PES, i.e., adiabatic 
and configuration-fixed (diabatic) constrained approaches. In the adiabatic approach, one 
always occupies the orbital that results in the lowest configuration energy. However, the 
PES obtained in adiabatic calculations may correspond to different configurations wher¬ 
ever there is orbital crossing. To obtain the PES for a given configuration, the concept 
of the so-called “parallel transport” should be used, which enables one to decompose the 
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ground state PES into separate parts uniquely characterized by the quantum numbers 


of the occupied orbits 


21 


22| . This conhguration-hxed constrained appro ach is often 


referred to as the diabatic constrained approach in the literature 


21 


22 


23|. 


The diabatic constrained approach is useful to investigate level crossings which lead 
to flip-flop situations. It is also useful in the Generator-Coordinate-Method (GCM) cal¬ 
culations in cases that excited conhgurations are included. Although the so-obtained 
excited states cannot always be identihed as physical states since the solutions of the cor¬ 
responding equations violate the variational principle, the diabatic constrained approach 
may serve as a basis to qualitatively understand the excited states and the interplay be¬ 
tween the excited states and the ground state. To describe these more quantitatively, one 
certainly should go beyond the mean-held level. 

Conventionally, the diabatic constrained method has been combined with the non- 
relativistic approaches, such as shell model, Skyrme or Gogny Hartree-Fock or Hartree- 
Fock-Bogoliubov methods, in understanding the properties of atomic nuclei, including 
high-spin states 23(], level crossings 2^, fusion and hssion processes 25(], etc. For exam¬ 
ple, the interesting phenomenon of shape coexistence 2^ can be studied using PES in 
complementary to the conventional interpretation of many particle many hole excitations 


27|. In the RMF model, the diabatic constrained method has not received enough at¬ 


tention compared with that of the adiabatic constrained method and the non-relativistic 
ones. Therefore, the purpose of this paper is to develop the diabatic constrained approach 
within the RMF model and the doubly magic nucleus ^‘^^Pb is chosen as an example for 
demonstration. 

This paper is organized as follows: Section II contains an outline of the RMF model. 
The adiabatic and diabatic constrained PES in ^°®Pb, together with the diabatic single 
particle spectra of the spherical ground state, are analyzed in Section III. The results are 
briefly summarized in Section IV. 


II. THEORETICAL FRAMEWORK 

The relativistic mean held model describes the nucleus as a system of Dirac nucleons 
which interact in a relativistic covariant manner via the meson helds. The meson helds 


3 














considered in most applications of the RMF model are the scalar-isoscalar sigma meson, 
vector-isoscalar omega meson, vector-isovector rho meson and the photon: the sigma me¬ 
son provides the long-range attractive interaction, the omega meson provides the short 
and mid range repnlsive interaction, the rho meson is responsible for the isospin depen¬ 


dence of the nuclear force while the photon accounts for the electromagnetic interaction. 


The Lagrangian density used in our present calculation is of the following form: 



where ^|J is the Dirac spinor and "0 = m, m^, and nip are the nucleon, a, uj 

and p meson masses respectively, while ga, g 2 , 5 ' 3 , guj, C 3 , gp, and ^ jAn = 1/137 are the 
corresponding coupling constants for the mesons and the photon. The held tensors of the 


vector mesons (a; and p) and of the electromagnetic helds take the following form: 


(2a) 

(2b) 

(2c) 


Ffc ^ g|.^x _ g^xl". 


With the classical variational principle, one can obtain the coupled equations of motion, 
i.e., the Dirac equation for the nucleons and the Klein-Gordon type equations for the 
mesons and the photon. In the most general case, these equations are very difhcult to 
solve; therefore various symmetries must be utilized to simplify them. 

In the case of an axially symmetric system, the projection of angular momentum on 
the z-axis, D, and the parity, vr, are good quantum numbers. Here, we consider only even- 
even nuclei, therefore time reversal symmetry is conserved. To solve the RMF equations 
for an axially symmetric system, we expand the Dirac spinors and the Boson helds with 
the eigenfunctions of a deformed harmonic oscillator potential in cylindrical coordinates 








^ ig^{z, ^ 


( 3 ) 


4 






20 oscillator shells have been used to expand both the Fermion helds and Boson helds to 
guarantee convergence. 

The potential energy surface is obtained through the constrained calculation. More 
specihcally, the binding energy at certain deformation is obtained by constraining the 
quadruple moment (Q 2 ) to a given value /i, i.e., 

{H') = (H)+ \c(m - p.)\ (4) 


with 


{Q 2 ) — {Q2)n + {Q2)p, 


( 5 ) 


where {Q 2 )n,p = (2r^P2(cos6*))„^p = {2z^ — — y‘^)n,p and C is the curvature constant 

parameter. The more often used deformation parameter (32 is related with the expectation 
- - 3 

value {Q 2 ) by (Q 2 ) = _ Ar'^P 2 , r = {Rq = 1.2 fm) and A is the mass number. 

vStt 

By varying fi, the binding energy at different deformation can be obtained. In principle, 
one has to follow a multidimensional energy surface. However, for the sake of efficiency, 
a one-parameter line, e.g. 02, is often used in microscopic constrained RMF calculations 


29j and the small effects of other deformations, such as the 


due to the self-consistency 
hexadecupole deformation. 

In order to obtain the diabatic PES with a hxed conhguration, the occupied orbits are 
determined by the so-called “parallel transport”, i.e.. 


('/'i(«)l'/’i(9+ Ag))|A,^ci « 


( 6 ) 


where i and j enumerate all the single-particle states of two adjacent conhgurations. In 
such a way, the original conhguration at q can be tracked and the corresponding PES can 
be obtained as a function of the deformation 301. 

In principle, if Aq is small enough, the conhguration at g -|- Aq should be the same 
as that at q. However, due to numerical difficulties, condition in Eq. ([6]) can not be 
rigorously implemented. Therefore, the following procedures are adopted in our study: 
First, the wavefunction and occupation number of every state at the initial conhguration 
q are recorded. Second, for each state i at conhguration q, we search all the states of 
conhguration q + Aq for state j that has the largest overlap with state i. It is to be noted 
that state i and state j have the same quantum numbers Once state j is determined. 
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the occupation number of state i will be transferred to state j. This procedure is repeated 
till the occupation number of each state at configuration q + Ag fixed. In case that many 
single particle states with the same symmetry are very close to each other, we adjust the 
step size of Ag to avoid the mismatch of the single-particle conhguration. 


III. RESULTS AND DISCUSSION 


In the preceding section, we described how to implement the diabatic constrained 
method within the relativistic mean held model. In this section, we would like to compare 
the results obtained with this method and those with the usual adiabatic constrained 
one. For this purpose, we study the potential energy surface of ^°®Pb. This nucleus is 
chosen due to the following considerations: First, its properties can be well described 
by the RMF model; second, as a doubly magic nucleus, the pairing correlation can be 
safely ignored. The latter is particularly important because the pairing correlation can 
complicate the issue dramatically 2^. Taking into account the pairing correlation, such 
as using a constant gap pairing scheme, does not qualitatively alter our results and the 
corresponding conclusions. However, in the presence of the pairing correlation, due to 
the fractioned occupation number, it is not convenient to track conhguration changes and 
discuss the relevant physics. Therefore, only the results without the pairing correlation 
are presented. It should be mentioned that the parameter set used is PKl, which has 
been carefully adjusted to reproduce both the nuclear matter saturation properties and 
the ground-state properties of hnite nuclei with a microscopic recipe for the center of mass 
correction 31|. 


In Fig. [T], the PES of ^osp]^ obtained from adiabatic (open circle) and diabatic (solid 
line) calculations are plotted as functions of the deformation parameter (32 ■ Surprisingly, 
it is found that even for ^°^Pb, which has a rigid spherical shape at ground state, the 
adiabatic PES is somewhat complicated, which can be easily seen from the many broken 
regions on the PES. The fact that there exists un-converged region on the adiabatic 
PES has been known for a long time. Over the years, it has been argued that it might 
originate fronr (i) the abrupt change of configuration at a certain point due to the no 
crossing rule [291]; (h) change of the mean helds due to varied conhguration [30(|; and/or 
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(iii) the increased mixing of all low-lying states due to their near degeneracy 3^ . On the 
other hand, the diabatic method is known to be able to connect broken regions existing 
in the adiabatic method. This can be clearly seen in Fig. [TJ In addition, from Fig. [1], as 
expected, one can see that the diabatic results (solid lines) can reproduce the adiabatic 
results (open circles) very well if they both exist at the same region. 

It should be noted that several local minima which do not exist on the adiabatic PES 
appear on the diabatic one. These minima, which could be viewed as excited states, are 
denoted as stars and labelled by capital letters in Fig. [1] Since in the diabatic method, 
the single-particle levels of each conhguration can be unambiguously determined (see 
the preceding section), we tabulate the corresponding configuration of each new local 
minimum in Table [Tl Each minimum is labelled by a bracket {Ex,f32), where is the 
excited energy relative to the ground state A and (32 is the corresponding deformation 
parameter. The conhguration of each state is given with respect to its core listed in the 
second row of the table. The conhguration of state F relative to that of state E is not 
given simply because the diherence between these two conhgurations is too large to ht 
into the table. We would like to point out a few more interesting things: 


• From Table [H it is seen that for most cases both neutrons and protons are excited 
as conhguration changes from one to the neighboring excited state. This might be 
related to the strong interaction between neutrons and protons. In this sense, a 
self-consistent calculation, such as the present one, is important to take care of this 
interaction and yields reliable results. 

• Due to the fact that diherent local minima have diherent conhguration, even the 
barrier between them is very small, they can still be well dehned, such as B and A, 
F and G 21]. 

In Fig.l, the diabatic potential (solid curve) in Fig. 1 show discontinuity around 
the deformations (3 ~ 0.4, which is due to the absence of the other diabatic curves 
with diherent conhgurations. If more careful diabatic constrained calculations with these 
conhgurations were done, the diabatic curves will become continuous and the two values 
at some quadrupole deformations, e.g., (3 ~ 0.4, would belong to one or two of these 
diabatic curves. 
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To provide a microscopic explanation of Fig. [H we plot the single particle spectra 
of near Fermi level obtained from the diabatic calculation in Fig. [2l Solid (dot- 

dashed) lines represent positive (negative) parity, filled (open) circles label the occupied 
(unoccupied) states, and the original configuration is that of the spherical case A (see 
Fig. [I]). The level crossing regions are highlighted by rectangles. From the left panel 
of Fig. [2], one can easily see those high j levels originating from 1 / 19/2 come down while 
those low j levels, such as 3 si /2 and 2 ^ 3 / 2 , go up as the nucleus becomes more oblate. 
At around P 2 ~ —0.1, the first level crossing happens between 3 si /2 and one lhg /2 state. 
The same is true for neutrons, in particular, the first level crossing happens between 
3 pi /2 and one I/ 11/2 state at \P 2 \ slightly larger than that where the first proton level 
crossing occurs. It shows that a diabatic calculation is important to obtain the smooth 
evolution of the single-particle levels as a function of deformation. Otherwise irregularity 
of the single-particle levels will appear due to the occupation of the lower orbits in the 
adiabatic calculation. We should point out that here proton and neutron level crossings 
occur almost simultaneously mainly due to the fact that their corresponding Fermi levels 
are close to each other in the case of ^°^Pb. It may not be the case for either proton 
or neutron rich nuclei where the proton and neutron Fermi surfaces are quite different. 
Beyond the crossing points, high j levels continue coming down while low j levels continue 
going up, which results in a very steep potential energy surface. This is the case for most 
diabatic PES’s displayed in Fig. [TJ 

Finally, we plot the adiabatic single particle levels in Fig. [3l They look similar to 
those shown in Fig. [2] but differ in two aspects. First, the diabatic single particle levels 
are continuously distributed as functions of (32 but the adiabatic single particle levels are 
not. Their broken pattern is the same as that of the adiabatic PES shown in Fig. [T] The 
underlying reason is, as we have mentioned above, that the sudden jump of the adiabatic 
PES is mostly due to the abrupt change of the corresponding configuration [^. For 
example, the crossing of 3 pi /2 and liii /2 at (32 ~ —0.1 results in the broken PES shown in 
Fig. [H Second, the diabatic single particle levels are occupied according to the dictated 
conhguration but the adiabatic single particle levels are so occupied to give the largest 
binding energy. 


IV. SUMMARY 


A diabatic constrained relativistic mean field approach is proposed to calculate the 
potential energy surface of the nucleus and applied to Although both adiabatic 

and diabatic constrained calculations yield almost the same ground state PES, the dia¬ 
batic one has the advantages that it enables one (i) to decompose the segmented PES 
obtained in usual adiabatic approaches into separate parts characterized uniquely by dif¬ 
ferent configurations; and (ii) to define the single particle orbits at each deformation by 
their quantum numbers unambiguously determined from their counterparts at the spheri¬ 
cal configuration. Thus, physics behind the adiabatic PES can be understood more clearly 
in the diabatic picture. 

This work is partly supported by the National Natural Science Foundation of China 
under Grant No. 10435010, 10221003, and 10447101, the Doctoral Program Foundation 
from the Ministry of Education in China. 
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TABLE I: Configurations of the five local minima (B, C, D, E, G) on the diabatic potential 
energy surfaces of ^ggg pjg, [J]), Each minimum is labelled by a bracket {Ex, (52) where 

Ex is the excited energy relative to the ground state A and P 2 , the corresponding deformation 
parameter. The configuration of each state is relative to the core as listed in the second column. 


State 

Core 

Configuration 

TT 

u 

B (7.63,-0.11) 

A 

(3'Sl/2) ^( 1 ^ 9 / 2 )^ 

(3pi/2) ^( 1 * 11 / 2 )^ 

C (15.63,-0.28) 

B 

( 2 ^ 3 / 2 ) ^( 1 * 1372 )^( 2 / 7 / 2 )^ 

( 3 ^ 3 / 2 ) ^( 259 / 2 )^( 1 / 15 / 2 )^ 

D (21.19,-0.39) 

C 


( 2 / 5 / 2 ) ^(2h9/2)^(l/l3/2)^ 

E (21.13,0.21) 

A 

( 3 S 1 / 2 ) ^( 2 ^ 3 / 2 ) ^( 1 ^ 9 / 2 )^ 

( 3 ^ 1 / 2 ) ^( 2 / 5 / 2 ) ^( 259 / 2 )^( 1 * 11 / 2 )^ 

G (19.19,0.69) 

F 

( 157 / 2 ) ^( 1 ^ 11 / 2 ) ^( 1 ^ 9 / 2 )^( 1 * 1372 )^ 

( 3 ^ 3 / 2 ) ^( 2 / 7 / 2 ) ^ 1 * 11 / 2 )^ 1 / 15 / 2 )^ 





FIG. 1: (Color online) Potential energy surfaces of ^'^^Pb obtained in adiabatic (open circles) 
and diabatic calculations (solid lines). The local minima on the corresponding diabatic PES are 
denoted as stars and labelled by different capital letters, whose configuration, binding energy 
and deformation are tabulated in Table. HI 
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FIG. 2: (Color online) Proton (left panel) and neutron (right panel) single particle states of ^osp]^ 
obtained in diabatic calculations as functions of the deformation parameter (32- Solid line and 
dot-dashed line represent positive and negative parity states, respectively. Filled (open) circle 
denotes whether the corresponding state is occupied (unoccupied). The original configuration 
is that of the ground state. For the sake of clarity, only those states near the fermi surface are 
shown; and the level crossing points are highlighted by rectangles. 



FIG. 3: (Color online) Proton (left panel) and neutron (right panel) single particle states of 
208pb obtained in adiabatic calculations as functions of the deformation parameter (32 ■ 
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